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As an alternative to the paradigmatic fragmentation problem of a single object crushed into a 
great number of pieces, we survey a large collection of identical bodies, each one randomly split into 
two fragments only. While some key features of usual fragmentation are preserved, this minimal 
approach allows for closed analytical results on both, shape abundances and mass distributions for 
the fragments, with robust power-law regimes. All the results are compared to numerical simulations. 
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The generic designation complex system fits few sub- 
jects as well as fragmentation. Though the interest in this 
field grew from specific examples [JJ , it now spreads from 
basic physics to material science and geology, including 
multifragmentation of nuclei Q, fracture propagation in 
diverse materials Q, and earthquakes [4j. The typi- 
cal fragmentation problem consists in the catastrophic 
crushing |5j] of a body due to the sudden appearance of 
fracturing stresses, the main goal being to determine the 
mass distribution of the resulting myriads of pieces. The 
outcome of such a process strongly depends on a large 
number of factors and their mutual interactions. Among 
them, the mechanicalproperties of the material, e. g., 
brittle or plastic Q, its effective dimensionality 
the mechanism of fragmentation (explosion, impact by a 
projectile, uniform compression, etc) , and the magnitude 
of the energy input [10]. This multitude of variables gen- 
erates a rich pool of data which has often been observed 
to follow power laws. A fully analytical approach to this 
class of problems is far beyond our present capabilities, 
both, due to mathematical complexity and to the lack 
of a complete knowledge on the physical responses that 
lead to crack nucleation and propagation under different 
conditions. 

Although the reduction of an intricate problem to a 
simpler, but non-trivial one, is a common modus operandi 
of physics, this was not fully carried on in the field of frag- 
mentation (see however the classical works by Gilvarry 
[HI]). Instead of considering one object broken into an 
enormous number of pieces we focus on an set of identi- 
cal bodies minimally segmented, i. e., each one randomly 
split into two fragments only. Suppose you take a tile and 
throw it from a modest height so that the scenario where 
a single crack develops is the most likely. By repeating 
this procedure with a large number of elements and se- 
lecting those minimally fragmented, one can obtain sta- 
tistical information on the shape and mass distribution of 
the resulting shards in the same way as it is done in usual 
fragmentation problems. In fact, the experimental setup 
presented in Ref. [12| on fiat bodies laterally impacted, 
seem appropriate to generate this kind of ensemble in a 
controlled way (for spherical bodies see [HI). A further 



FIG. 1: (color online) Minimally segmented tiles. In (a) the 
fragments are a triangle and a pentagon, while in (b) two 
trapezoids are obtained. In (c) we represent a minimal frag- 
mentation. The variables I and <j> axe shown along with the 
adopted labeling of the sides. 



example of a spontaneous process that gives rise to min- 
imal fragmentation is the propagation of a fracture on a 
tiled floor. 

We will be concerned with 2D objects, namely, rectan- 
gles with side lengths a and 6, defining an aspect ratio 
7 = b/a. We notice that minimal segmentation is natu- 
rally associated to a low energy process, the cracks are 
neither branched nor too irregular [see figs. [Ha) and 
QJb)]. For this reason, in our geometric statistical model, 
they are represented by a random segment for each ele- 
mentary rectangle. We must make the meaning of these 
segments completely unambiguous. As discussed in sev- 
eral works, there is no a priori geometric statistical pre- 
scription to define the partitioning of a body with di- 
mension greater than one. Different rules lead to quite 
distinct results, a fact that is strikingly illustrated by 
the Bertrand Paradox Since we are interested in 
situations where the cracks either originate or cross the 
boundaries of rectangles without any bias and isotropi- 
cally, our basic assumptions are: the segmentation line 
has a uniform probability to cross any interval in the 
boundary of the rectangle, the angular distribution of 
the segments being also uniform. Therefore, in order to 
generate a line we define two random variables I and <f> 
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uniformly distributed in the intervals [0, 2a + 2b] = [0, L] 
and [— 7r/2, tt/2], respectively. For definiteness, <j> is given 
as the angle between the segment and the normal direc- 
tion to the side selected by the variable I [see Fig. He)]. 
This construction encompasses the desired physical in- 
puts and, at the same time, employs the minimum pos- 
sible set of variables. 

It has been recently suggested that geometric proper- 
ties, in opposition to mass distributions, may be a more 
stringent test to theoretical models [Hj]. In this regard 
the shape of the resulting fragments can be classified 
according to the way a single line can split a rectan- 
gle. The two alternatives are: (i) the intersection points 
with the border are on neighboring sides or (ii) the seg- 
ment cross opposite sides. Within our approximation, 
scenario (i) produces as fragments a triangle plus a pen- 
tagon, while situation (ii) leads to two trapezoids, all 
irregular in general. We begin by addressing the follow- 
ing question: what is the relative frequency of events (i) 
and (ii) as a function of the aspect ratio 7 for a large 
number of events? We denote all (statistically) relevant 
crossings as follows: the segment jk intersects sides j 
and k, with j, k = 1,2,3,4, [Fig. [He)] the associated 
probabilities given by Pjk("f)- Scenario (i) occurs if the 
segment is of the type 13 or 24. For situation (ii) the 
possibilities are 12, 23, 34, 41. Thus, the probability 
of obtaining a crack that connects neighboring sides is 
mapped into P n = P12 + P23 + P34 + Pax- From the 
symmetry of the problem it is evident that these prob- 
abilities are equal, so P n — 4Pi2. To calculate P12 we 
note that 12 may arise in two ways. First, the point 
drawn is on side 1 (i £ [0, a]) and the angle satisfy 
cot -1 [6/ (a — I)] < 4> < tt/2, which gives the angular in- 
terval A0 = tt/2 - cot _1 [&/(a - I)] = tan" 1 [6/(a - I)]. 
Since the angular distribution is uniform, the probability 
of obtaining 12 in this first case is 

F( 7 ) = i Ttan- 1 f_L 
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The second possibility of getting 12 is I £ [a, a + b] and 
— tt/2 < (f> < cot -1 [a/ (I — a)], that yields the same result 
as in (JTJ, with a ■<-> b, i. e., 7 —> 1/7, as expected. We 
then get P12 = F( 7 ) + F(l/j), leading to 

Pn = 4F( 7 ) + 4F(l/ 7 ) . (2) 

The probability of two opposite sides traversed by the 
random segment (two trapezoids) is P = 1 — P n . In Fig. 
[5] we show these probabilities as functions of 7. The con- 
tinuous curves represent P n and P a . The open triangles 
and squares correspond to numerical experiments with 
10 6 segmentation events (for each point). Note that P n 
attains its maximum value for 7 = 1, where 

Pn = \ + — * 0.72 . (3) 



This probability remains larger than P a for 7 < 3.83, the 
aspect ratio of the rectangle highlighted in Fig. [2j 




FIG. 2: (color online) Probabilities of obtaining two trape- 
zoids (P ) and a triangle plus a pentagon (Pn) as functions of 
7. For 7 w 3.83 these probabilities are equal. 

Note, however, that these results say nothing about 
the size of the fragments. Obtaining an experimental 
plot compatible to Fig. [5] is not as simple as it may 
appear. In order to get the right frequencies, tiny frag- 
ments (which can be easily overlooked) must be counted 
as much as larger ones. The point is that the density of 
very small triangles diverges as the fragment mass ap- 
proaches zero, a signature of usual multifragmentation. 
This might not be expected at a first glance, since our 
model seems to produce mainly sizable fragments. This 
naturally leads us to the central fragmentation problem 
of mass distributions. In experiments it is much easier to 
record the accumulated probability V> (m) for a fragment 
to have a mass greater than m. However, for convenience 
we will use the complementary accumulated probability 
V(m) = 1-V > {m): 



V(m) 



dm'p(m') 



(4) 



where p(m) is the associated probability distribution. Let 
2M = ab be the total mass of the rectangle (without loss 
of generality we set the superficial density of mass to 1). 
Since in each event we have a minimal segmentation, for 
every fragment of mass m there is another one of mass 
2M-to, that is, V(m) = V>(2M-m) = l-V(2M-m). 
In particular, for m = M we must have V(M) = 0.5. 
Given this symmetry, we only need to consider values of 
m between and M, meaning that, in each event, it is 
sufficient to collect the smaller of the two pieces, which 
is either a triangle or a trapezoid, to get all information. 

We proceed by defining the auxiliary quantity V(l,m), 
the accumulated probability for a fixed value of I, i. e., 



V(m) 



V(l,m)dl 



(5) 



To determine V(l, m) we must specify the mass m' < m G 
[0,M] of an arbitrary smaller fragment in terms of the 
random variables I and <f>. Let us suppose once more that 
side 1 is selected, I £ [0, a]. In this case, if the fragment 
is a triangle, then ml = (1/2) (a — Z) 2 cot</>, while for a 
trapezoid we get m! — b(a — I) — (l/2)& 2 tan0. Note 
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also that if m < b(a — l)/2, only triangles are possible. 
Given m and I the limiting angle for triangular fragments 
is 4>iim — tan -1 [(a — I) 2 /2m}. Again, since the angular 
distribution is uniform, the probability is proportional to 
the allowed angular interval A(f> = n/2 — (j>iim- Thus, we 
get 



for multifragmentation of planar objects 14j . This be- 
havior is completely dominated by triangular fragments 
since the probability of getting a trapezoid becomes van- 
ishingly small for ji — > 0. The final regime is observed 
for 7 7^ 1, becoming more evident for greater values of 7. 
For massive fragments we get 



V(l,m) = - 

7T 



7T , / (a — I) 2 

- - tan -1 v n ' 
2 V 2m 



(6) 



for < m < b(a - l)/2, that is, < I < a - 2m/b. 
Accordingly, if b(a — I)/ 2 < m < ab/2 we must have a 
trapezoid, leading to 



V'{l,m) = - 
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From these results we obtain the following accumulated 
probability with I restricted to side 1: 



V(l,m)dl + - 
L 



V'(l,m)dl 



(8) 



The full accumulated probability is then V(m) — 
Ti{m) +T 2 (m) +T 3 {m) +V A {m) = 2V X (m) + 2V 2 (m) = 
2Vi(m) + 2T'i 4r ^ b (m), where the equalities are justified 
by direct symmetry arguments. The final result can be 
compactly written as 
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The transition from PT|) to (fT2")) is due to a change in ef- 
fective dimensionality. To understand this, consider the 
minimal fragmentation of a thin, uniform, rod with unit 
size and mass. Let / be a random variable uniformly 
distributed on [0, 1], the value of I defining the point at 
which the rod breaks. It is evident that the mass distri- 
bution of smaller fragments directly follows the distribu- 
tion of I itself, i. e., p(l) oc 1° p(fi) cx The fact is 
that, especially for higher aspect ratios, only small frag- 
ments realize the 2D character of the elongated rectangle, 
for larger fragments the behavior of the probability den- 
sity is virtually the same as in the ID case. See Fig. 
[3l where the ratio m/M is nearly the same in both de- 
picted situations. For a fragment to be considered large 
its characteristic size should overstep the length of the 
shorter side of the rectangle, i.e., its area must exceed a 2 
in comparison to the total area ab (or a 2 /2 in comparison 
to M), which leads to fi ~ a 2 /ab = 1/7, so, for /1 > I/7, 

All these features are 



,(9) 



p is essentially given by ([T2|) . 



where Co = /i/7, (i = M7j an d M = m/M is the dimen- 
sionless mass with /x 6 [0, 1]. The normalization is such 
that 7^(1,7) = 1. Inverting Eq. (g]) and using Eq. © we 
get the following probability distribution 
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Therefore, we obtain all information on the ensemble in 
a fully analytical way. We find the density p to be es- 
sentially a composition of two power laws connected by 
a stretched dimensional crossover. First, in the limit of 
very small fragments, which we call the dust regime, we 
find that p asymptotically diverges as 
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(11) 



showing that the model leads, chiefly, to small fragments. 
Interestingly enough, this same exponent has been re- 
ported in numerical simulations of a number of models 



FIG. 3: (color online) For an elongated body the ratio be- 
tween the fragment mass and total mass is nearly the same as 
in the ID case. The square delimited by the dashed line gives 
the minimum area that characterizes a "large" fragment. 
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FIG. 4: (color online) Log-log plots forp(/i). The continuous 
curves refer to Eq. (|10l) and the open circles to numerical 
simulations. The power-law regimes are evident in the four 
panels concerning the depicted rectangles (7 = 1,3,5, 10). 



4 



shown in the log-log plots of Fig. [4] where the density 
p is plotted against fi for aspect ratios of 1 (a), 3 (b), 
5 (c), and 10 (d). The continuous curve represents the 
analytical result of Eq. (|10p and the open circles refer to 
numerical samplings with 5 x 10 6 minimal fragmentation 
events. The dust regime is shown to behave according 
to the power law given by Eq. (fTTj) . while for large frag- 
ments a flat region, Eq. (fT2|) . is observed for /i > 1/7 
(dashed vertical lines). 

To verify whether these exponents are robust with re- 
spect to variations in the model, especially, those which 
would make it more realistic, let us consider two relevant 
points. First, the triangle in Fig. (Ha) is smaller than 
the actual fragment. Assuming that the contrary situa- 
tion is equally likely, we repeated our numerical calcula- 
tions setting the new mass as m new = m(l + 5), where 
S is a uniform random variable over a symmetric inter- 
val [— e, e]. The results (not shown) are indistinguishable 
from those of the original model. Second, although it 
has been often used in the literature, the isotropy condi- 
tion may not be strictly followed, for example, the crack 
may have a higher probability to propagate in the vicin- 
ity of the dashed line in Fig. [He). We then assumed 
that ef> obeys a parabolic probability distribution, given 
by (3 — K(f> 2 for <fr € [— 7r/2, tt/2], with the normalization 
condition irj3 = 1 + 7t 3 k/12, and found that the exponents 
1/2 and are robust, the crossover region being shorter. 
In Fig. [5]/3 and n are such that the probability to find <j) 
in the interval [tt/2, tt/2 — d<fi] is 50% smaller than that 
of finding <f> within [0, d<fi}. 
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FIG. 5: (color online) Same as in Fig. [4] for non-isotopic 
random segments. The continuous line represents Eq. (|10p 
without modification. 

In contrast, alternative definitions of the random vari- 
ables themselves lead to very distinct results, showing 
that our choice of a point on the perimeter is an essen- 
tial physical assumption. We made calculations with the 
choice of an interior point plus an angle and the results 
do not resemble any data on fragmentation. In particular 
no power-law divergence has been observed. 

In conclusion, we addressed a new class of fragmen- 
tation problems. At first, it might be thought that the 
typical power-law divergence for small masses would not 
show up. Contrary to this expectation, we obtained an 
exact exponent of a\ — 0.5 in the dust regime. A second 
power law has been observed for large fragments with 
a-2 = 0. These exponents are insensitive to the aspect 



ratio of the rectangles, and robust with respect to varia- 
tions in the probability distributions used to generate the 
fragments. In addition, we determined the shape abun- 
dances of the fragments for arbitrary aspect ratios. It is 
expected that this kind of minimal approach may stimu- 
late new experimental realizations and theoretical efforts. 
The presented model do not intend to describe all quan- 
titative features of multifragmentation, even because the 
exponents are scattered over a large range of experimen- 
tal values. However, the simple model presented here 
describes experimental facts like the absence of a char- 
acteristic size in the small mass regime and transitions 
due to effective dimensionality. To better describe the 
range of exponents of multifragmentation, one possibility 
is to consider hierarchical models based on generations of 
minimal segmentation events, the number of generations 
being related to the energy input. 
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